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In this paper, we study the flux of electrons and positrons injected by pulsars and by annihilating 
or decaying dark matter in the context of recent ATIC, PAMELA, Fermi, and HESS data. We 
review the flux from a single pulsar and derive the flux from a distribution of pulsars. We point 
out that the particle acceleration in the pulsar magnetosphere is insufficient to explain the observed 
excess of electrons and positrons with energy E ^ 1 TeV and one has to take into account an ad- 
ditional acceleration of electrons at the termination shock between the pulsar and its wind nebula. 
We show that at energies less than a few hundred GeV, the flux from a continuous distribution of 
pulsars provides a good approximation to the expected flux from pulsars in the Australia Telescope 
National Facility (ATNF) catalog. At higher energies, we demonstrate that the electron/positron 
flux measured at the Earth will be dominated by a few young nearby pulsars, and therefore the 
spectrum would contain bumplike features. We argue that the presence of such features at high en- 
ergies would strongly suggest a pulsar origin of the anomalous contribution to electron and positron 
fluxes. The absence of features either points to a dark matter origin or constrains pulsar models in 
such a way that the fluctuations are suppressed. Also we derive that the features can be partially 
smeared due to spatial variation of the energy losses during propagation. 
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I. INTRODUCTION 

The nature of dark matter (DM) remains one of the 
most interesting problems in cosmology and astrophysics. 
Recently, several cosmic ray experiments reported higher 
than expected fluxes of electrons and positrons at ener- 
gies between 10 GeV and 1 TeV 0, i, S, 0, Q • One inter- 
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also a viable possibility. The main purpose of this work 
is to identify an observational signature which might dif- 
ferentiate between these two models. 

Pulsars are known to produce and accelerate electrons 
and positrons in their magnetosphere However, as 
we will argue below, this acceleration is insufficient to 
explain the observed excess. Moreover, the spectrum of 
particles in the magnetosphere is further modified at the 
termination shock between the pulsar and the interstellar 
medium (ISM) or the supernova remnant (SNR) which 
may contain the pulsar wind nebula (P WN) . It is impor- 
tant to emphasize that termination shocks are observed 
around middle-aged pulsars (e.g., Geminga 
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J1747-2958 [23|) and not just young pulsars like the Crab. 
Most young energetic pulsars are also surrounded by a 
PWN (for a review see, e.g., 2J|) which is powered by 
the continuous emission from the pulsar and remains in- 
side the shell or envelope produced by the initial explo- 
sion (a discussion of various regions surrounding a pulsar 
can be found in, e.g., 25|). Inside a PWN, the electrons 
and positrons are confined by the nebula's magnetic field 
for a long period of time before escaping into the ISM. 
Since a pulsar loses the vast majority of the spin-down 
energy while its PWN still exists, we assume that most 
of the electrons and positrons injected by pulsars spend a 
significant amount of time inside a PWN before reaching 
the ISM. This has two main consequences: 

• The spectrum of electrons and positrons injected 
by a pulsar into the surrounding ISM is not the 
spectrum of particles inside the magnetosphere (as 



assumed by many authors, e.g |15l.ll7l.l26|). but the 
spectrum of particles that escape the PWN into the 
surrounding ISM. In this paper, we assume that 
this is the same as the electron/positron spectrum 
inside the PWN when it is disrupted, which we es- 
timate using the observed broadband spectrum of 
these objects. 



Since the hfetime of a PWN {t < 100 kyr; [2||) 



is generally much smaller than the typical propa- 
gation time {t > 100 kyr), the electrons observed 
at the Earth come from PWNe that no longer ex- 
ist, whereas electrons inside existing PWNe cannot 
reach us. Since the variability of PWNe properties 
is very large, we cannot predict the electron flux 
from a pulsar that has already lost its PWN, even 
if we fully know its current properties (e.g age, po- 
sition, spin-down luminosity). 

Therefore, the best one can do is to use the currently 
observed PWNe to derive a statistical distribution of 
their properties. This distribution can be used to find 
either the average fiux of electrons expected from pulsars 
or, by assigning the PWNe random properties according 
to the distribution, the typical electron flux observed on 
Earth from all pulsars. 

The observed spectrum on Earth of electrons and 
positrons injected by pulsars is also strongly dependent 
on propagation effects. In particular, the observed cutoff 
in the flux of electrons from a pulsar can be much smaller 
than the injection cutoff due to energy losses ("cooling") 
during propagation. We define the cooling break, i?br(i)j 
as the maximal energy electrons can have after propa- 
gating for time t. Since - as stated above - the typical 
electron propagation time is much larger than the life- 
time of a PWN, we can assume that a pulsar is a delta- 
function source Q ~ d{t) and the propagation time for 
electrons from this pulsar can be estimated by the pul- 
sar's age t. If the cooling break is at a lower energy than 
the injection cutoff from a PWN, £'br(i) < £'inj, then 
the observed break is the cooling break which depends 
on the age of the pulsar but is independent of the in- 
jection cutoff. Since the cooling break is much steeper 
than the injection break, the existence of several pul- 
sars with i?br ^ ^'inj sufficiently close to the Earth such 
that the propagation time of electrons they inject into 
the ISM is less than their age will result in a sequence of 
steps or bumps in the spectrum. At high energies, only 
a few pulsars will satisfy this criterion, so these steps are 
expected to be well separated and therefore observable. 
At lower energies, we expect many pulsars to contribute, 
causing these steps to be averaged together and resulting 
in a smooth spectrum. We argue that the presence of 
significant steps, or bumps, in the electron spectrum at 
high energies would strongly suggest this excess is gen- 
erated by pulsars, since the fiux from the dark matter 
is expected to be smooth with a single cutoff. If these 
features are not observed, a pulsar explanation of the ob- 
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served excess is possible if there are no young energetic from pulsars and dark matter 



pulsars in the vicinity of the Earth with ^ i?inj 27 1 
or there are considerable spatial variations in the energy 
losses of these particles as they diffuse through the ISM. 
The amplitude of these fluctuations also depends on the 
relative contributions of the backgrounds and the pul- 
sars. If backgrounds dominate at high energies, these 
fluctuations may be undetectable. 

This paper is organized as follows: In Sec. |lTl we review 
the propagation of electrons in the ISM. We estimate the 
typical propagation time and distance for electrons and 
positrons of a given energy and show that injection of 
electrons and positrons by a pulsar can be approximated 
by delta-functions in space and time. Using the proper- 
ties of pulsars and their PWN derived in Appendix[Al we 
derive the temporal evolution of the observed flux from 
a single pulsar at a given distance. In Sec. IIIIl we study 
the flux from a distribution of pulsars, calculating the 
average expected flux and comparing this result to the 
estimate flux from pulsars in the Australia Telescope Na- 
tional Facility (ATNF) catalog Q - demonstrating that 
at energies E < 300 GeV this flux is well approximated 
by the average curve while at higher energies there are 
significant deviations and bumplike features due to the 
cooling breaks. We then compare both spectra with the 
ATIC, Fermi-LAT, and PAMELA data. In Sec. Ml we 
review the flux from dark matter. In Sec. |Vl we present 
our conclusions. We argue that the flux from dark matter 
is likely indistinguishable from the flux from a single pul- 
sar or from a continuous distribution of pulsars. Thus, 
unless there are additional features at high energies that 
point to the contribution from several pulsars, it may 
be impossible to tell whether the excess is due to dark 
matter or pulsars. 

The paper also contains a few appendixes: In Ap- 
pendix [X] we give a general review of pulsars and pulsar 
wind nebulae. In Appendix [B] we derive the smearing 
of the cooling breaks due to spatial variability of en- 
ergy losses. In Appendix [C] we study some constraints 
that current electron and positron data put on the pul- 
sar models. 



II. SINGLE PULSAR FLUX 

The flux of cosmic ray electrons at the Earth depends 
on both the spectrum of injected electrons and the prop- 
erties of the ISM. First, we review the propagation of 
electrons in the ISM, and then derive the expected flux 



A. Properties of the interstellar medium 

The ISM contains a magnetic held with a strength on 
the order of 3/iG [2^. Since the corresponding Larmor 
radius for a 1 TeV electron is small, ''"l = f§ < 
pc, electrons are expected to mostly follow the ISM's 
magnetic field lines. Because the ISM magnetic field has 
random fiuctuations, electrons propagate along on a ran- 
dom path. The corresponding diffusion coefficient for 
relativistic particles is 30 1 

5 

(1) 



DiE) ^D,[§- 



where (5 = 0.3 - 0.6 and Do = (3 - 5) x lO^^cm^ s^^ for 
Eq — \ GeV. [The typical mean free path for a 1 TeV 
electron rt ^ D{E)/c > 1 pc is much larger than the 
Larmor radius tl < 10^"^ pc] We find it convenient to 
express the diffusion coefficient in terms of the energy 
rather than the magnetic rigidity, R = p/q^ where p is 
the momentum and q is the charge of the particle. In our 
case, the two definitions are equivalent. 

As electrons propagate in the ISM, they lose energy. 
For electrons with energy E > 5 GeV, the dominant loss 
mechanisms are synchrotron radiation and inverse Comp- 
ton scattering off cosmic microwave background, infrared 
(IR), and starlight photons. The corresponding energy 
losses are 



E = ^b{E) = -boE^ 



(2) 



where 60 = 1-6 x 10 ^''GeV ^s ^ for the local density of 



photons 3l| and B = 3 fiG. Before we present a formal 
solution to the propagation equations, we define a few 
characteristic numbers. In estimations, it is convenient to 
represent the parameters in the units of pc = 3 x lO^^cm 
and kyr = 3 x 10^'^s. In this paper we will usually use 



bo = 1.6 X 10 
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GeV^'s- 
2^-1 



5 X 10-^GeV^^kyr"^ and 



Do = 3 X lO^^cm^s-i = lOOpc^kyr"^ with S = 0.4. 

By integrating Eq. ([2]), we find the energy loss in terms 
of the electron travel time is 



1 1 

= bnt. 

El Eo " ' 



(3) 



where Eq is the initial energy of the electron and Ei is 
the energy at time t. The cooling break, defined as the 
maximal energy an electron can have after traveling for 
time t, is 

1 
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(4) 
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The characteristic travel time is therefore 

t > 100 kyr for E <2 TeV. (5) 

The characteristic distance an electron travels before 
cooling to energy E is the diffusion distance x'^^^f — 
4D{E)t, where t=^, 



a;diff<5kpc for i;>10GeV. 



(6) 



B. Green function for diffusion-loss propagation 

In general, the evolution of the energy density p of elec- 
trons moving in random paths and losing energy can be 
described by the following diffusion- loss equation 2^ 32 1 



0(x,i?,<), (7) 



where Q = dN/{dEdt(fx) is the energy density of elec- 
trons injected by the source. In principle, one can also 
take into account reacceleration, convection, and decays 
(collisions), but for electrons with i? > 10 GeV these 
contributions can be ignored. 

The general solution to Eq. ^ is found in (j^lfist. 
To solve Eq. ([7]) for a general source, one introduces the 
Green function G(x, E, t; xq, Eq, to) which satisfies 

(8) 



d{^~^o)SiE-Eo)S{t-to). 
Then, the solution to ([7]) is 
p(x, E, t) = 

j d^xo j dEo J dto G(x, E, t; xq, Eq, h) 
•g(xo, Eo, to). (9) 
The Green function can be derived as follows. One can 



define the variables t' ~ t — t and A 
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where 



T EE t{E, Eq 

A(S, Eo) = 



-Eo 



dE' 



E b{E'y 
^" D{E')dE' 



E 



b{E') 



(10) 
(11) 



The variable t' {t, E) is invariant with respect to the dif- 
ferential operator dt — h{E)dE- In fact, D~^{E){dt — 
b{E)dE) = d\ and Eq. ^ becomes the usual diffusion 
equation in A and x. The Green function is then [s^jsS 

1 1 (x-xq)^ 

G'(x, E, t; xo, £^0, h) = 



b(E) (47rA)3/2 



Equation ([7]) and the above Green function have a few 
limitations. Both the ISM magnetic field and density of 
IR and starlight photons vary in space. Consequently the 
diffusion coefficient and the energy loss function depend 
on the coordinates: D = D{E, x), b = b{E, x), and there 
is no simple analytic solution to Eq. ([7]). In Appendix [Bl 
we calculate corrections to the predicted e'^e~ spectrum 
at the Earth due to spatial variations in the energy loss 
function. 

C. Flux from a single pulsar 

With the general Green function in hand, one can find 
the density of electrons at any point in space for any 
source. In this Section, we derive the expected flux of 
electrons and positrons produced by a single pulsar. 

The distances to pulsars are sufficiently large that we 
can assume that pulsars are point sources. We also as- 
sume that most of the pulsars' rotational energy is lost 
via magnetic dipole radiation 3J] which eventually trans- 
forms into the energy of electrons and positrons 

-2 



gpui.ar(x,s,i) = Q(£;)i ( 1 



t 



9{t)d{x), (13) 



where t is the pulsar age and x is its position. 9{t) is the 
step function that ensures Qpuisar = for t < 0. Note 
that the pulsar spin-down time scale r in this formula 
and the variable introduced in (fTO|) are unrelated. We 
review the derivation of this formula in Appendix El 

At late times (t ^ r) , the spin-down luminosity scales 
as t~^. Consequently, most of the energy is emitted dur- 
ing t ^ T. The pulsar spin-down time scale t < 10 kyr is 
much smaller than the typical electron propagation time 
t > 100 kyr. Consequently, we can take the limit r ^ 0, 
which results in 

-2 



S{t). 



(14) 



For pulsars with a significant PWN, the time dependence 
in Eq. (|13p does not describe the escape of electron and 
positrons from the PWN into the ISM. However, for most 
pulsars the lifetime of the PWN {t < 20 kyr; j2^) is 
much smaller than the typical propagation time i > 100 
kyr [Eq. Thus, the detailed time dependence of the 
escape is not significant and the delta-function approxi- 
mation in Eq. ^4]) is still valid. 

We assume that the energy spectrum of particles in- 
jected into the ISM Q{E) by a given pulsar has the form 



■S{t-to^T)eiEo-E). 



(12) 



(15) 
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where n is the injection index and M is the injection 
cutoff. We denote the initial rotational energy of the 
pulsar by Wq and define rj to be the fraction of this energy 
deposited in the ISM as e+e~. The total energy emitted 
in e'^e~ is then 



J Q{E)EdE = rjWo. 



For n < 2, this gives 



r(2 - n)M- 



2-n ■ 



(16) 



(17) 



The index n and the cutoff M of the electron spectrum 
can be derived from the broadband spectrum of a PWN, 
and may vary significantly between pulsars. We argue in 
Appendix [X] that reasonable values for energetic pulsars 
are n = 1.5 ± 0.5 and M - 100 GeV - 10 TeV. 

The overall normalization is more difhcult to derive 
because neither the initial rotational energy nor the con- 
version efficiency are known for most pulsars. Currently, 
models can only estimate VFo to an order of magnitude 
with significant theoretical uncertainties 35|. In Ap- 
pendix 1^ we derive that, assuming a constant pulsar 
time scale r = 1 kyr for all pulsars in the ATNF cata- 
log [3] ) the distribution of pulsar initial rotational ener- 
gies Wo = 1 OP erg satisfies a log-normal distribution with 
p « 49 and ap k, 1 which gives the average Wq ~ 10^° erg. 
If we use r = 10 kyr, the same analysis gives p ~ 48, 
CTj, ~ 1, and Wq w 10'*^ erg. If the age of a pulsar is 
known independently, then the initial rotational energy 
can be estimated more robustly. For example, the Crab 
pulsar is associated with the SN1054 supernova explosion 
and has r « 0.7 kyr and Wq ~ 5.3 x 10*^ erg. 

Let us now discuss the conversion coefficient rj. The 
energy density near the surface of the pulsar is domi- 
nated by the magnetic field and the spin-down luminos- 
ity is dominated by the magnetic dipole radiation. In 
most PWNe the energy density is believed to be parti- 
cle dominated, i.e. at large distances from the neutron 
star most of the energy outflow has been converted to 
particles and 77 '--^ 1 at this stage (see, e.g., [i^ for a dis- 
cussion of the Crab PWN). However, these particles do 
not immediately escape to the ISM but are trapped in- 
side the PWN by its magnetic field where they can lose 
a significant fraction of their energy. In Appendix \^ 
we estimate 77 ~ 0.1 due to cooling of the particles be- 
fore they escape into the ISM. Based on the discussion 
above, we find that the average energy in electrons and 
positrons rjWo ^ 10^^ erg is reasonable. This value is 



model dependent and can vary greatly from one pulsar 
to another. 

The density of electrons propagated from a pulsar to 
the Earth can be found by substituting the source func- 
tion Q(x, E, t) into Eq. (O 



p(x, E, t) 



b{Eo) 



1 



h{E) (47rA)'V2 



QiEo), 



(18) 



where parameter A is defined in Eq. (|lip and Eq is the 
initial energy of the electrons that cool down to E in time 
t 

E 



Eq 



1 - Ebot 



(19) 



The density in Eq. has a cutoff at the cooling break, 

= jij, since Q{Eq) for Eq — > 00. 

For a density p of relativistic particles, the flux is de- 
fined as 



F = —p. 



(20) 



The time evolution of the flux from a single pulsar is 
shown in Fig. [TJ At early times, the electrons have not 
had enough time to diffuse to the observer and the flux is 
exponentially suppressed. At later times, the flux grows 
until the diffusion distance AD{E)t is similar to the 
distance from the pulsar to the Earth. After that the 
flux decreases as the electrons diffuse over a larger vol- 
ume. The cutoff moves to lower energies due to cooling 
of electrons. 

For energies much smaller than the cooling break, we 
can neglect energy losses. In this case, Eq E, \ k, 
D{E)t and Eq. ^ reduces to 

1 



p(x, E, t) = 



_e TomTQ{E). (21) 



{4TrD{E)t)y 

Assuming that < AD{E)t, the flux for i? < is 



F{E) 



At: {ATiDQtf/^ 



E 



-n-4<5 



(22) 



In general, the flux that we add to the backgrounds to fit 
the data can be effectively parametrized by three num- 
bers: the normalization, the index at low energies, and 
the cutoff energy. From the right hand side of Eqs. (fT7|) 
and we find that these three parameters correspond 
to at least 8 parameters describing the pulsar and the 
ISM. In particular, the propagated index Ua is a linear 
combination of n and 5, Ua — n + |(5, the propagated 
cutoff is i^cut = and the normalization depends on 77, 
Wo, M, Do J and t. In order to fix this degeneracy, as a 
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FIG. 1: Time evolution of e"'"e~ flux on the Earth from a 
pulsar at a distance of 1 kpc with rjWo = 3 x 10*® erg, an 
injection index n — 1.6, and an injection cutoff M = 10 TeV. 
The diffusion and energy losses are described in Sec. Ill Al 
We assume the delta-function approximation for the emission 
from the pulsar, Q{x,E,t) = Q{E)5{x)5{t). The flux from a 
young pulsar (the 3 kyr curve on the right) has an exponential 
suppression because the electrons have not had enough time 
to diffuse from the pulsar to the Earth. The cutoff moves 
to the left due to cooling of electrons and becomes sharper. 
After reaching a maximal value, the flux decreases since the 
electrons diffuse over a large volume. 



matter of convenience, we choose Dq = 3 x lO'^^cm^s ^, 
S = 0.4, bo = 1.6 X 10-i^GeV"^s-i, Wo = 10^° erg, and 
M = 10 TeV. With this choice, our fit to the e+e~ data 
wiU determine n, 77, and t. If some of the parameters 
are known independently, e.g., the propagation model, 
the energy losses, the age of the pulsar etc., this ap- 
proach becomes more constrained and more predictive. 
As shown in Fig. [21 the expected flux from a pulsar 
with rjWo w 3 X 10"*^ erg, n = 1.6, distance 0.3 kpc, and 
age 200 kyr reproduces the positron fraction measured 
by PAMELA and is a good fit to the cosmic-ray electron 
spectrum measured by ATIC, Fermi, and HESS below 
^ 1 TeV. This suggests that the anomaly in the e~''e~ 
flux could be due to a single pulsar. However, given the 
considerable number of known nearby, energetic pulsars 
[2^, it is unlikely that the flux from any single pulsar is 
significantly larger than the flux from all such pulsars. 
In the next section, we will derive the expected flux of 
electrons and positrons from a collection of pulsars. 



III. FLUX FROM A COLLECTION OF 
PULSARS 

In this section, we derive the e+e^ flux from a con- 
tinuous distribution of pulsars and compare it with the 
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FIG. 2: Electron and positron flux from a single pulsar to- 
gether with a primary background ~ E~^'^ and a secondary 
background ~ E~''''^ . The pulsar is at a distance of 0.3 kpc. 
It has ryVKo = 2.2 x 10*® erg, age of 200 kyr, and injection 
index and cutoff n — 1.6 and A/ — 10 TeV, respectively. The 
propagation parameters are described in Sec. Ill Al The cutoff 
Af 3> 1 TeV results in a significant bump around 1 TeV which 
is consistent with the ATIC data. For a smaller injection cut- 
off M ~ 1 TeV, the flux from the pulsar takes the form of a 
power law with an exponential cutoff that can be used to flt 
the Fermi and PAMELA data (see, e.g., [s^]). 



predicted fiux from the pulsars in the ATNF catalog [2 



A. Flux derivation 

We assume that pulsars are homogeneously distributed 
in the galactic plane and are born at a constant rate Nb 
35|. The "continuous" distribution of pulsars is defined 
as the average of all possible realizations of pulsar dis- 
tributions. This results in a source function constant in 
time, localized in the vertical direction, and homogeneous 
in the galactic plane 



Qdi«tr(x,S,i) = Jo E-'' e-^ Siz) 



(23) 
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FIG. 3: The expected spectrum from the continuous flux dis- 
tribution and that from pulsars in the ATNF catalog pulsar 
[2^ . The latter is calculated using 77 = 0.065, n = 1.5, and 
a pulsar time scale t = 1 kyr for each pulsar. This last fact, 
in conjunction with its spin-down age and current spin-down 
luminosity, is used to calculate each pulsar's initial rotational 
energy through Eq. (I29p . We also use the value of the propa- 
gation parameters given in Sec. Ill Al Several hundred pulsars 
contribute below 300 GeV and the continuous distribution 
provides a good approximation for these energies. Above 300 
GeV, there are only ^ 10 contributing pulsars, and the ob- 
served flux in this energy range is strongly dependent on their 
individual properties. The reason for the significant discrep- 
ancy between these two curves above 2 TeV has to do with 
the actual local distribution of pulsars versus the averaged 
flux seen by many observers in the Galaxy, as discussed in 

SedmH 
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FIG. 4: The predicted flux from pulsars in the ATNF cat- 
alog calculated using the same procedure as in Fig. |3] but 
accounting for spatial variations of energy losses as described 
in Appendix|B] The assumed backgrounds are the same as in 
Fig.E 



with the normalization constant 



Jo 



(24) 



r(2 - ji)7\f 2-" Ag^i ' 

where A^^\ is the area of the galactic plane. Since the dif- 
fusion distance of these electrons is significantly smaller 
than the distance from the Earth to the edge of the galac- 
tic plane [il] (xdiff < 10 kpc), we can neglect the effects 
of having an edge at a finite distance. 

Using the general Green function in Eq. (fT^ , the flux 
of electrons from this distribution is 

F = 

£- y" d^xo y" dEn J dta G(x, E, t; xo, -Bo, to) 
•Q(xo, So, to). (25) 
Integrating over to and x, we obtain 

/•OO -1 

C / ,„ 1 



FiE) 



dEo 



--JoEo^e- 



(26) 



where A is defined in Eq. pT|) . This flux can be rewritten 

as 

c Jo 



Fdistr{E) 



47r V47r6oA 



Ie_ £;-«-('5+i)/2^ (27) 



where 



Ie ^ 



dx 



1 



(28) 



for example, \i E <^ M, S — 0.4, and n = 1.5, then 
JjL « 3. 

M 

As in the case of a single pulsar flux, the number of 
parameters we need to fit the data is much smaller than 
the number of parameters characterizing the flux from 
a collection of pulsars. In this case, the index of the 
observed flux and the normalization can be found from 
Eq. (|27p. For example, the index of the flux at low en- 
ergies Ua ~ n + {1 + S)/2. Formally, the cutoff in this 
case is equal to the injection cutoff M, but for an actual 
distribution of pulsars the expected cutoff is lower and 
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is determined by the age of the youngest pulsar within 
the diffusion distance from the observer - as derived in 
Sec. lIII B1 If we break the degeneracy by picking a partic- 
ular propagation model, we can constrain the properties 
of the pulsar distribution. The opposite is also true: by 
choosing some properties of the pulsars one can constrain 
the properties of the ISM - as demonstrated in Appendix 

o 

In order to break the degeneracy we fix the ISM proper- 
ties as in Sec. Ill CI To calculate the flux from the pulsars 
in the ATNF catalog we use the following toy model. We 
assume that every pulsar has injection index n = 1.5 and 
conversion efficiency rj = 0.065 (these values are chosen 
to fit the low energy electron and positron data in Fig. 
S]). We choose an injection cutoff M = 10 TeV for ev- 
ery pulsar (for smaller values of M the features at high 
energies will be less sharp, since the injection cutoff is 
not as abrupt as the cooling break) . In order to estimate 
the initial rotational energy, we assume that for each pul- 
sar, the spin-down time scale is r = 1 kyr. Then we use 
Eq. (jA4[) to express the initial rotational energy £q = Wq 
in terms of the current spin-down luminosity £ and the 
pulsar age t: 

W^^E-. (29) 

T 

The result in shown in Fig. [31 and the relative normal- 
ization between this spectrum and that of the continuous 
distribution described above depends on the pulsar birth 
rate A^b, or, to be more precise, on the local value of the 
pulsar birth rate. In order to have a good agreement 
between the two curves for energies 30 - 300 GeV, we 
require A'b ~ 1-8 kyr~^, assuming a Milky Way radius 
Rgai — 20 kpc [35]. For energies below 30 GeV, we find 
that the main contribution to e^e~ flux comes from the 
pulsars with age t > 10 Myr. These pulsars typically 
have a very low spin-down luminosity and therefore are 
difficult to observe (in the ATNF catalog there are very 
few pulsars with the spin-down luminosities E < 10^^ 
erg). In Fig. [H we apply to this spectrum the Gaussian 
smearing expected to result from spatial variations in en- 
ergy losses depending on the path of the electrons - as 
derived in Appendix |B] As one can see, this provides a 
very good fit to the PAMELA, Fermi, and HESS data - 
but does not reproduce the ATIC bump. 

Determining the flux of e+e^ from the actual distribu- 
tion of pulsars using a more realistic model is extremely 
difflcult because every pulsar has its own independent 
parameters (e.g., Wq, rj, and r). Thus, we may choose 



several thousands of parameters in order to flt less than 
a hundred of data points (which can be fltted by a flux 
parametrized by three parameters only). Moreover, as we 
discussed in the Introduction, these thousands of param- 
eters refer to PWNe sufficiently old that their electrons 
have had enough time to diffuse to the Earth. These 
PWNe have already disappeared and therefore cannot 
be observed directly - making it impossible to directly 
constrain these parameters observationally. The large 
number of pulsars and the impossibility to derive their in- 
dividual properties suggest a statistical method is needed 
to study the e+e~ fiux they produce. At small energies, 
a lot of pulsars contribute to the observed flux on Earth 
and therefore the properties of an individual pulsar are 
unimportant. In this case, the flux should be well ap- 
proximated by some average curve - as demonstrated in 
Fig. [3l In this estimate we included all pulsars with ages 
t > 15 kyr and use the delta- function approximation of 
source functions, (5(x, E, t) = Q{E)S{-x.)S{t) (Sec. 
The choice of the lower cutoff on the age of the pulsars is 
motivated by the fact that young pulsars, such as the Vela 
pulsar, usually have a PWN and therefore their electrons 
have not escaped yet into the ISM. At high energies only 
a few young pulsars contribute and the deviation from 
the average curve may be large. The presence of features 
at high energies may serve as a signature of a collection 
of pulsars that can distinguish them from a dark matter 
or single pulsar origin for these electrons. 

B. Statistical cutoff 

As shown in Fig. [3l the expected flux from a contin- 
uous distribution of pulsar increases with energy until a 
break at ^ M — 10 TeV, whereas the predicted flux 
from pulsars in the ATNF database has a cutoff at 2 TeV. 
This discrepancy is due to the rare events when a young 
pulsar is very close to the observer. Since electrons lose 
energy during propagation, high energy electrons must 
come from young pulsars and the cutoff energy is deter- 
mined by the age of the youngest pulsar sufficiently close 
to the observer so that the electrons have enough time 
to diffuse through the ISM. We will call the average such 
cutoff a "statistical" cutoff. 

To estimate the statistical cutoff, we consider a collec- 
tion of pulsars and choose an observation point. The sta- 
tistical cutoff at this point is the maximal cooling break 
energy for the flux from these pulsars. In this distribu- 
tion, the youngest pulsar whose electrons can reach the 
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FIG. 5: Statistical cutoff as a function of the diffusion in- 
dex and the birth rate of pulsars in the Galaxy. The cutoff 
in e^e~ flux from pulsars is determined by the age of the 
youngest pulsar within the diffusion distance from the Earth. 
The average such cutoff is a universal quantity that depends 
on the properties of ISM (the energy losses and the diffusion 
coefficient) and on the pulsar birth rate, but it is insensitive to 
the properties of the injection spectrum from the pulsars. We 
assume Do = lOOpc^kyr"^ and 6o = 5 x IO~''GeV"^kyr"\ 

observation point has an age T and diffusion distance 
R. For a given pulsar birth rate iVb, we estimate Mstat 
by demanding that there is at least one pulsar within R 
younger than T. Therefore, we have a system of three 
equations for the three unknowns R, T and Mstat- 



1 

i?2 = 4Z?(M,tat)T, 

N,T— = 1. 
Solving this system of equations, we find 



V ^0 ^gal 



(30) 
(31) 
(32) 



(33) 



Assuming i?gai = 20 kpc, Dq = 10 ''kpc^ kyr ^, and 
bo = 5x 10"*^GeV"^kyr~\ we get 



(4 X lO^TVb) 



GeV, 



(34) 



where A^b is in units of kyr^^. In Fig. [Sj we show the 
statistical cutoff as a function of iVb and the diffusion 
index 6. This calculation should be viewed as a rough 
estimate, with the actual flux from the distribution of real 
pulsars having a cutoff that differs by as much as an order 
of magnitude. Additionally, it is possible that current 
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FIG. 6: The flux from a contiimous distribution of pulsars. 
The parameters are chosen to fit the Fermi and PAMELA 
data points, riWo~6-5 x 10** erg and n — 1.5. In this plot, 
instead of the injection cutoff AI = 10 TeV, we use the statis- 
tical cutoff A/stat = 1 TeV. The backgrounds are the same as 
in Fig. (2] The propagation parameters are described in Sec. 

EH 



data are missing a feature at high energies {E > 2 TeV) 
due to poor statistics. A comparison between the flux 
from a continuous distribution of pulsars with Mstat — 1 
TeV and the current data is shown in Fig. [6l 

We note that Eq. ([55]) can also be used to find the 
cutoff in the primary background if we assume that it is 
generated by the supernova explosions. For instance, for 
the supernova rate in the Milky Way A^sn = lOkyr^^ and 
6 = 0.4, it gives the cutoff in the primary background 
around 3 TeV. Using the same reasoning as above one 
may expect some features in the spectrum of the primary 
electrons at several TeV. Below ^ 1 TeV we do not expect 
significant fluctuations in the primary background and 
the presence of the features should be interpreted as the 
signature of pulsars. 
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The source function of e+e coming from annihilating 
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FIG. 7: Flux from DM model in with the annihilation 
chain x + X — > + fli — > 2e+ + 2e", Mdm = 1 TeV and a 
boost factor BF ~ 500. The primary background is ~ E~^ '' 
and the secondary background is ~ 

IV. FLUX FROM DARK MATTER 

In this section we briefly review the e~''e~ production 
from annihilating (decaying) DM and derive that, for a 
large class of DM models, the expected flux has the form 
of a power law with a universal index n — 2 at energies 
E <C A/dm. If we neglect gradients in the DM density 
near the Earth, then we approximate any DM contribu- 
tion as originating from a constant, homogeneous source, 
which from Eq. ^ gives 



PiE) 



1 



HE) 



Q{E')dE'. 



(35) 



This equation has an interesting property that for any 
Q{E) ~ E^ with A: > — 1, the integral is saturated at the 
upper limit, which in this case is the mass of the DM 
particle A/dm- For energies E <^ A/dm, we can neglect 
the dependence on E resulting from the lower limit of 
integration so the index of the electron flux is determined 
by the index of the energy loss function b{E) ^ E'^ . 



dark matter is 



38| 



QiE) = -nliav)-, 



(36) 



where n-^ is the dark matter number density, (av) is 
the thermally averaged annihilation cross-section, and 
dN / dE is the number density of electrons and positrons 
produced per annihilation event. Here we assume that 
the DM particle is its own antiparticle otherwise there is 
an extra factor of 1/2 in Eq. ([36|) . For this source func- 
tion, the flux of electrons and positrons from annihilating 
DM is 



F{E) = 



1 



Stt h{E) 



n\{av) 



M 



dN 

Ie 



{E') dE'. 



(37) 



If the integral in this equation is saturated at E^ < A/dm, 
then for E <^ E^, the integral is insensitive to the changes 
of the lower integration limit and can be approximated 
by a constant 



L± = 



f^^ dN 



(38) 



where /g± is the average number of electrons and 
positrons produced in an annihilation event. In this case, 
the only energy dependence in F{E) is from b{E), so 
F{E) ~ E~^. The discussion of the universality of index 
n = 2 with respect to the choice of DM models and DM 



halo profiles is further discussed in [39| (see also [40|[41| 
for an earlier discussion of the effects of DM substruc- 
ture). 

An important difference between the DM and pulsar 
models is that the dark matter flux in Eq. ((37|) has signif- 
icantly fewer free parameters than the corresponding flux 
from pulsars. In fact, if we assume that the energy losses 
in the ISM are well understood and the energy density 
of dark matter is fixed from the cosmological considera- 
tions, there are only two free parameters. A/dm and (cru), 
with the specific DM model providing /e± . For a given 
DM model. A/dm is then fixed by the cutoff energy in 
the observed spectrum and the cross section (av) is fixed 
by the normalization of the flux. The index of the flux 
is not parametrically independent, n fa 2. This index is 
insensitive to the choice of DM model or the DM pro- 
file in the host halo but may change significantly in the 
presence of a large DM subhalo 39|. As an example, 
we use the DM model in [ij^ with the annihilation chain 
X + X — ' (t' + 4' — * 2e+ + 2er. DM with the current 
estimated energy density of — 0.3 GeVcm""^ requires 
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(crw)o = 3.0 X 10 ^^cm^s ^ at freeze out. One can as- 
sume that the current cross section is larger by a boost 
factor (BF). To fit the ATIC and PAMELA data, we set 
A/dm = 1 TeV, which requires a BF ^ 500 to reproduce 
the observed normalization (Fig. [7]). 

For a decaying DM model, Eq. ([57]) would be replaced 

by 



FiE) 



1 



47r b{E) Td Je dE' 



M 



dN 



-dE', 



(39) 



where is the life-time of the DM particle and / ^ dE 
is the number of electrons and positrons produced per 
decay. If we take the same number density and the mass 
of DM particles as above, then 



-^5x 10-27s-i. 

Td 



(40) 



These estimates agree with the analysis of [4^ [43| j44| . 



V. CONCLUSIONS 

In this work, we analyzed the flux of electrons and 
positrons from a single pulsar, from a continuous distri- 
bution of pulsars, from pulsars in the ATNF catalog and 
from dark matter. Depending on the model parameters 
and pulsar properties, they all can adequately fit either 
the Fermi and PAMELA data or the ATIC and PAMELA 
data. One of the most important question is whether it 
is possible to distinguish among these possibilities. 

In Fig. [8] we compare the expected e+e~ flux from a 
single pulsar (Sec. |TT|, pulsars in the ATNF catalog (Sec. 
mil) , a continuous distribution of pulsars (Sec. IIII[) . and 
DM (Sec. lIVp . We have chosen the parameters of the 
models such that the fluxes have the same value at 100 
GeV, similar indices at low energies, and a cutoff at 2 
TeV. At energies below ^ 300 GeV, the fluxes are very 
similar. We also do not expect to see any differences be- 
tween these models in the positron ratio below 300 GeV, 
the upper limit for charge identification in PAMELA. 

Above 300 GeV, there are substantial differences 
among the e'^e^ spectrum predicted for these models. 
However, it should be noted that the sharpness of the 
cutoff for a single pulsar and for DM is strongly model 
dependent. If the injection cutoff for a pulsar is ^ 1 
TeV, then the cutoff in the observed flux from a single 
pulsar can be much smoother than if the injection cut- 
off was higher than the cooling break, in which case its 
spectrum is indistinguishable from that predicted for a 
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FIG. 8: The fluxes from annihilating dark matter, from a sin- 
gle pulsar, and from a continuous distribution of pulsars can 
be made similar, depending on the parameters of the models. 
The flux from a collection of pulsars may have significant devi- 
ations from a continuous curve. This property can be used to 
distinguish the pulsars from the sources producing a feature- 
less spectrum. The flux from pulsars in the ATNF catalog is 
the same as in Fig. U The single pulsar has the age t = 100 
kyr, distance 0.3 kpc, rjWo = 9.2 x 10** erg, and n = 1.6. 
The continuous pulsar distribution has A'^b = 1.8 kyr~^, 
riWo = 6.5 X 10*** erg, n = 1.5, and A/^tat = 2 Tev. The dark 
matter model is the same as in Sec. IIVI but with A^dm = 2 
TeV and BF = 2000. The ISM properties are the same as in 
Sec. UlAl 



continuous distribution of pulsars. For the DM flux we 
show a model with only one intermediate particle in the 
annihilation-decay process. If there are more steps in the 
annihilation-decay process, then the flux has a broader 
cutoff and, again, may be impossible to distinguish be- 
tween either a single pulsar or continuous pulsar distri- 
bution origin. Thus, given the significant uncertainties 
in the pulsar and DM models, it is unlikely that better 
observations alone can distinguish between a single pul- 
sar and dark matter origin of anomalous e"'"e~ flux 26| 
(a similar conclusion was obtained in [3] 45|). 

The flux from a discrete collection of pulsars does have 
a few distinctive features at high energies. The height of 
these features is model dependent and may be within the 
error bars of current observations. The presence of these 
features requires the existence of a few young, nearby, en- 
ergetic pulsars with an injection cutoff ^ 1 TeV. Conse- 
quently the absence of such features in the observed e~'"e~ 
spectrum could mean that all young pulsars whose elec- 
trons have reached the Earth have had cutoffs < 1 TeV - 
a strong constraint on the properties of PWNe since, as 
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we discuss in Appendix^l the PWN around the Vela pul- 
sar has a cutoff in the electron and positron spectrum at 
an energy ^ 1 TeV. An additional smearing of the bumps 
can be due to spatial variations of the energy losses and 
the diffusion coefficient. Our general conclusion is that 
the current electron and positron data are not sufficient 
to distinguish between the pulsars and the dark matter 
and that independent measurements, e.g., the spectrum 
and morph olog y of the diffuse galactic gamma-ray back- 



ground 45|, l4a | , may be necessary in order to decisively 



distinguish a pulsar and DM origin of e+e excess. 
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APPENDIX A: REVIEW OF PULSARS 

In this appendix we review the emission of electrons 
from pulsars. We assume that this emission is powered 
by the pulsar's loss of rotational energy. Pulsars are be- 
lieved to be rotating neutron stars with a strong surface 
magnetic field Q, and magnetic dipole radiation is be- 
lieved to provide a good description for its observed loss 
of rotational energy. A pulsar loses its rotational energy 
on a characteristic decay time r defined as 



(Al) 



where £q and £q are the initial rotational energy and the 
initial spin-down luminosity, respectively, which in the 
magnetic dipole radiation model are equal to 



£o — -/Oq, 



So = 



B^R^ sin^ a 
6^ 



"0' 



where flo is the initial angular velocity, R is the radius of 
the pulsar, B is the strength of the surface dipole mag- 
netic field, and a is the angle between the rotation axis 
and the magnetic field axis. 

If the energy loss is due to magnetic dipole radiation, 
then 



mn = - 



B^R^ sin^a 
6c^ 



Integrating the energy loss equation we get 



£(t) = £o-[l 



t 



(A2) 

(A3) 
(A4) 



As a result, the pulsar angular velocity satisfies 

-(t + r). (A5) 



2h 



In a more general approach, the time evolution of the 
angular velocity is described as 



(A6) 



where k is the breaking index, which can be found by 
measuring the current Q, il, and 
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In this case, 



Q{t) = no 



T 



(A8) 



The magnetic dipole radiation corresponds to fc = 3. 

As an example, let us calculate the initial rotational 
energy of the Crab pulsar using the magnetic dipole ap- 
proximation and a general braking index. The Crab pul- 
sar is believed to have been produced during SN 1054 
supernova explosion. Consequently, the age of the pulsar 
is known exactly, t = 955 yr. In the magnetic dipole ap- 
proximation, we can use Eq. (jA5P and the current values 
of n and (l 3, 47 1 to calculate that its pulsar time scale 
T Ri 0.3 kyr. Assuming a mass 1.4 M©, radius R = 12 
km, and moment of inertia / — 1.4 x 10^^ g cm^ [j^ ]. 
we derive an initial rotational energy Wq « 3 x 10^° erg. 
Taking into account the measured value of O [3], the 
braking index of the Crab pulsar is k = 2.5, which gives 
T w 0.7 kyr and Wq w 5.3 x 10*9 erg. 

If the age of the pulsar is not known independently, it 
is impossible to determine r and Wq using its observed 
properties (from 17, 17, and f2 one can only calculate t+r). 
In the following, we estimate Wo by assuming t = 1 kyr 
for all pulsars. If so, the initial energy Wq = 8q can 
be found from Eq. (jA4[) by using the current spin-down 
luminosity £ 





£q = £t{1 + t/r) 



(A9) 
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FIG. 9: The left panel has 41 pulsars within 4 kpc from the 
Earth and younger than 300 kyr. Assuming the braking in- 
dex 3 and the pulsar time scale r = 1 kyr we find the average 
power p = Logio(Wo/erg) = 49.1 ± 1.1 with x jdoj = 1.2. 
For the right panel we select the pulsars from the ATNF cat- 
alog that have braking index 2 < fc < 10 (there are 20 such 
pulsars). Assuming the pulsar time scale r = 1 kyr, these 
pulsars have p = 48.9 ± 1.1 with jdoj = 0.3. For the pul- 
sar time scale r = 10 kyr, the left (right) selection of pulsars 
would have p = 48.2 ± 1.1(48.1 ± 1.2). Thus, for r = 1 kyr 
(10 kyr) the average initial rotational energy is Wo ~ 10^" erg 
(10« erg). 



This result is strongly dependent on the chosen value of t. 
For example, if r = 10 kyr, then an analogous calculation 
gives p « 48, dp w 1 and Wq « lO^'^erg. The estimations 
above agree with the analysis of |33 . It is worth noting 
that T likely varies between pulsars. 

To estimate ry, it is important to understand how the 
pulsar's magnetic dipole radiation is transferred to the 



For a general braking index /c, this formula takes the form kinetic energy of particles. Since £ decays as t"^, most 



=fr(l +</t)^ 



(AlO) 



Using this method to estimate Wq for all pulsars in the 
ATNF catalog results in the distribution shown in Fig. 
ini using both the magnetic dipole approximation (Eq. 
(IA9p ) and a general braking index method (Eq. (jAlOp ). 
In the magnetic dipole case we took all pulsars within 
4 kpc from the Earth and younger than 300 kyr. The 
reason is that older and more distant pulsars are less 
luminous and may not be observed for small spin-down 
luminosities, i.e., this introduces a bias towards more en- 
ergetic pulsars and shifts the distribution toward larger 
average Wq. For the general braking index, we used all 
pulsars with 2 < fc < 10. In both cases, Wq can be 
described by a log-normal distribution with the average 
•p — Logig(Wo/erg) w 49 and the standard deviation 
CTp w 1. The average initial rotational energy in this case 



Wq 



27r(T„ 



lO^e 



dp « 10^^ erg. (All) 



of the rotational energy is lost at early times. A young 
pulsar is surroimded by several layers [2^[23|. Nearest 
to the neutron star is the magnetosphere, which ends at 
the light cylinder i?Lc = c/H. The rotating magnetic 
field creates a strong electric field capable of both pro- 
ducing pairs of particles and accelerating them to rela- 
tivistic energies. These particles stream away from the 
light cylinder as a coherent "wind" that ends with a ter- 
mination shock separating the wind zone from the PWN 
which consists of magnetic fields and particles moving in 
random directions. The PWN in turn is surrounded by 
an SNR. A significant PWN exists only at the early times 
(t < 100 kyr 0). 

The spectrum of electrons and positrons in the magne- 
tosphere can be estimated using the spectrum of pulsed 
7-ray emission from a pulsar. For the Crab pulsar, model 
fits to the observed photon spectrum suggest that the 
spectrum of e+e~ pairs in its magnetosphere is well de- 
scribed by a broken power low with an index of 2.0 below 
ii^br ~ 2 GeV and an index of 2.8 between i^br and an 
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upper cutoff around 100 GeV [48[ (see also [49|). Parti- 
cles with this spectrum cannot reproduce the e~^e~ spec- 
trum observed on Earth, since a break in the injection 
spectrum at 2 GeV is too low to explain the ATIC and 
PAMELA results, and an index of 2.8 above 2 GeV also 
does not fit the data. Additionally, the pulsed emission 
from a pulsar only reflects the energy spectrum of the 
emitting particles in the emission region, which is not 
necessarily representative of the spectrum of particles 
that escape the pulsar magnetosphere along the open 
field lines and are eventually deposited in the ISM. 

These particles are further accelerated before they en- 
ter the PWN, most likely at the termination shock be- 
tween the magnetosphere and the PWN (for a review 



see, e.g., |50|). Once deposited in the PWN, they are 
trapped by the PWN's magnetic field until it is dis- 
rupted. Observationally, the spectrum of the electrons 
inside the PWN is found by analyzing their broadband 
spectrum, which at low photon energies (< 1 GeV) is 
dominated by synchrotron emission and at higher ener- 
gies (> 100 GeV) dominated by inverse Compton scat- 
tering of electrons off background photons [36]. From 
the radio spectrum of these objects, it is possible to con- 
strain the spectral shape of the low energy (GeV) elec- 
trons which dominate by number the electron population 
of a PWN. An average index of the electron and positron 
spectrum can be found using data from the publicly avail- 
able Catalogue of galactic SNRs [si] ]. where F-type (or 
"filled-center" ) SNRs are PWNe, S-type are the super- 
nova shells, and C-type SNRs are a combination of the 
two. There are 7 F-type SNRs with an average electron 
index np = 2a + 1 = 1.3 ± 0.3. For 21 C-type SNRs the 
average index is uq = 1.8 ± 0.4, while 168 S-type SNRs 
have ns = 2.0 ±0.3. It is clear that the spectrum of elec- 
trons in supernova shells is much softer (decreases faster 
with the energy) than the spectrum in PWNe. In order 
to explain the PAMELA positron ratio we need either F- 
type or, possibly, G-type SNRs because the S-type SNRs 
are produced by the initial supernova explosion and do 
not contain a significant number of positrons. 

The broadband spectrum of most PWNe shows a break 
between the radio and X-ray regimes, believed to corre- 
spond to a break in the electron and positron spectrum, 
most likely the result of synchrotron cooling. Gonverting 
the frequency of this break to an electron/positron energy 
requires knowing the strength of the PWN's magnetic 
field. An independent estimate of the magnetic field is 
available for those PWNe with detected inverse Compton 



emission, since this depends solely on the energy spec- 
trum of electrons and positrons in the PWN and known 
properties of the various background photon fields (e.g., 
Cosmic Microwave Background and starlight). The best 
studied example is the Crab Nebula (e.g., [Sa]), whose 
broadband photon spectrum suggests an electron spec- 
trum well described by a broken power law with an index 
n = 1.5 below Eq ^ 200 GeV and n = 2.4 between Eo 
and an upper cutoff Ecut ~ 10'^ TeV: the magnetic field in 
this PWN has a strength of i? « 2 x 10"'' G, resulting in 
a ratio of magnetic energy flux to particle energy flux of 
a < 0.01 36|. For PWNe whose broadband spectrum is 



not as well determined, the break energy in the electron 
spectrum is typically derived using the minimum energy 
assumption (a = 0.75 [54I). Using this method and the 



observational data provided in [52|, we estimate a break 
energy of ~ 3 — 300 GeV for the PWNe listed in this pa- 
per. It is important to emphasize that this procedure al- 
most certainly overestimates the magnetic field strength 
inside a PWN since, for most PWNe, a is believed to be 
(T ^ 0.75. In this procedure, the inferred break energy 



where i?pwn is the strength of PWN's 



Eb is Eb CX Spwn^, 

magnetic field. As a result, the true break energy of elec- 
trons and positrons inside the PWNe analyzed above is 
likely to be at least an order of magnitude higher than 
the derived value. 

The break energy in the electron/positron spectrum of 
a PWN is expected to vary considerably during the life- 
time of a PWN due largely to changes in the strength 



of the PWN's magnetic field (e.g., |53l, [SJ). There- 
fore, the break energy in the spectrum of electrons and 
positrons injected by the PWN into the surrounding ISM 
depends strongly on the evolutionary phase of the PWN 
when this occurs. During the initial free-expansion phase 
of the PWN's evolution (the Crab Nebula is the proto- 
typical example of such a PWN; [2J[), the break energy 
is expected to increase as ~ t^-^ [53|, [54 1. This phase 
of the PWN's evolution ends when it collides with the 
SNR's reverse shock, typically on the order of ~ 10^ yr 
after the supernova explosion. If this holds for the Crab 
Nebula, its current age of ^ 1000 yr and break energy 
of ~ 200 GeV suggests that, at the time of this colli- 
sion, the break energy will have risen to '--^ 8 TeV. The 
evolution of the break energy after this collision is more 
complicated and depends strongly on the properties of 
the central neutron star, progenitor supernova, and sur- 
rounding ISM (see [5^ for a more detailed discussion). 
There is observational evidence that the break energy of 
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older PWNe (> 10^ yr old) is considerably higher than 
that of the Crab Nebula and other young PWNe (~ 10^ 
yr). The most convincing example comes from a recent 
analysis of the broadband spectrum (radio to TeV 7 rays) 
of the Vela PWN (often referred to as "Vela X"), which 
suggests a break in the electron spectrum of ~ 67 TeV 
[sSj]. For the purpose of the work presented here, the 
exact value of the break is not important as far as it is 
bigger than « 1 TeV. 

Observations indicate that most PWNe are particle 
dominated: i.e., almost 100% of the spin-down luminos- 
ity is transformed into the energy of the particles after 
the termination shock. However, not all of this energy is 
eventually deposited in the ISM. We estimate this frac- 
tion, ?7, by first assuming that the spectrum after the 
termination shock is a power law Q{E) ^ E'"" with an 
index n < 2 and a cutoff ^ 10^ TeV. Then, the total 
energy in electrons is 



E-"EdE - E, 



2~n 



(A12) 



If, when the PWN is disrupted, the energy spectrum of 
electrons in the PWN is below the break at E'br and 
E~'^^ with Tib > 2 above the break, the total energy in 
electrons is saturated at E\^t- with 



Wfin ~ E, 



The efRciency is therefore 



br 



Ehi 

'eT 



(A13) 



(A14) 



For n = 1.5, E^, = 10 TeV and E^ = 10^ TeV, this 
gives the suggested = 0.1. This derivation should be 
viewed as an order of magnitude estimation. A more re- 
alistic calculation is extremely complicated and involves 
the knowledge of the PWN evolution and the actual spec- 
tra of particles inside a PWN. Recent work in this field 
does support an efficiency of 77 ~ 0.1 (e.g., Q)- 

It should be stressed that, apart from theoretical un- 
certainties, the parameters of the injection spectrum can 
vary significantly between pulsars. The initial rotational 
energy can differ by several orders of magnitude, while 
the index of the electron spectrum n can vary from 1 to 
2. In some cases, it is observed to vary inside the PWN 
of a single pulsar. The upper cutoff M, as we have seen 
in the examples of Crab and Vela pulsars, can vary at 
least between ~ 100 GeV and ~ 10 TeV. 



APPENDIX B: SPATIAL VARIATION IN 
ENERGY LOSSES 

As we have discussed, an important signature of the 
flux from the pulsars is the presence of a number of 
bumps at the cooling break energies 



bti 



(Bl) 



where ti's are the ages of the pulsars. The existence of 
these bumps is based on the assumption that the energy 
losses depend only on the travel time and not their path. 
In reality, the energy loss coefficient depends on the posi- 
tion, since the densities of the star light and IR photons 
vary in space. In this case, there is no simple solution 
for Eq. ([7]), though one can still find the average energy 
loss and its standard deviation by averaging the energy 
losses over random paths. 

As a useful simplification we will consider separately 
diffusion in space and energy losses. Our motivation is 
that a particle detected with energy E has an energy close 
to E during most of the propagation time (i.e. the cool- 
ing time from Eq to E is saturated by the final energy E). 
Consequently, the diffusion coefficient for all particles de- 
tected with energy E can be approximated by D{E). If 
so, the probability to propagate from a source at {xq, to) 
to an observer at [xi, ti) is given by the Green function 



G{xi, ti; xq, to) 



1 



[AttD{E)MY/^ 



g 4D(£;)At 



(B2) 



In order to find the energy loss averaged over paths, it is 
useful to rewrite this Green function in terms of the path 
integral 



G(xi, ti] xo, to) = j Dx{t) 



-S[x{t)] 



where the action is 

S[x{t)] 



AD{E) 



-x^dt 



(B3) 



(B4) 



with the boundary conditions x{ta) ~ x^ and x{ti) — xi. 

In general, the average of a functional ©[^(t)] over 
paths is 

_ !Dx{t)0[x(t)] e-^[-Wl 

- ^Dx{t) e-s[-it)\ ■ ^^^^ 

Integrating the energy loss 
dE 



dt 



-h{x)E' 



(B6) 
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along a path x{t), we find 
1 1 

The functional that we will study is 



0[x{t)] = / b{x)dt. 

J to 

The expression in the numerator of (|B5P is 



(B7) 



(B8) 



-S[x{t)] 



I Dx{t) I dt' b{x{t')) 

J J to 

dt' f Dx{t) b{x{t')) e-^[^(*)l. (B9) 



Num = 



If we define x' — x{t'), then all the paths can be repre- 
sented as a path from xq to x' , the integral over all x' 
and the path from x' to xi: 



Num = / dt' Dx{t) I dx' I Dx{t) 

J to ^ Xq J ^ x' 

■b{x{t')) e-S[-(*)] 

dt' I dx' G{xi, ti; x' , t') 



(BIO) 



■b{x') G{x', t'; xo, to). 



(Bll) 



The resulting expression resembles the first order pertur- 
bation theory: there is a propagation from xq to x' , an 
insertion of an operator at x' and a propagation from x' 
to Xi. 

The average energy loss can be estimated with (IBlip 
and the Green function in (jB2p . Taking i?o oo in (|B7p , 
we find that the average cooling break energy for a given 
pulsar is 

(-^) = {0) = { dt' f dx' G{xi, ti; x', t') b{x') 

-C/br J to J 

•G(x', t'- xo, to) )/G(a;i, t,; xo, to). (B12) 
The standard deviation is 

ao = ^/W)~W- (B13) 

The average (O^) for the functional (jB8|) can be com- 
puted analogously to (jB12p 

(©2) ^ 2G-\xi- xo) 

j ' dt' j dx' j ' dt" j dx" G{xi; x") 
■b{x") G{x"; x') b{x') G{x'; xq), 



where we assume that x' is at t' and x" is at t". The 
factor of 2 is the usual n! for the time-ordered path inte- 
grals. 

The relative standard deviation of the cooling break 
energy is 



AE _ ao 



(B14) 



Using the energy densities of starlight and IR photons 
from [si] ] we find the relative smearing in the energy 



AE 



E \ 



^1/3 



0.053- I I 

ITeV J 



(B15) 



At 1 TeV the smearing is about 5%, at 100 GeV it is 11%, 
and at 10 GeV it is 24%. The flux from the ATNF pulsars 
with this smearing is shown in Fig. H) At low energies the 
flux becomes very smooth but at high energies the bumps 
are still visible. We also notice that at high energies the 
relative width of the bumps is larger than the ratio in 
(jBlsp . Thus, even if the experimental energy resolution 
is about 10% - 15%, we should be able to see the bumps. 



APPENDIX C: CONSTRAINING PULSARS AND 
ISM PROPERTIES 

If we assume that the anomalies in Fermi, ATIC and 
PAMELA data are due to pulsars, we can use these re- 
sults to constrain the properties of ISM and pulsars. The 
problem is that the flux depends on both the properties 
of the ISM and the injection spectrum from pulsars. As 
we discuss in Sees. HIl and Hill the ISM can be described 
by three parameters Do, S, and bo and the injection from 
a distribution of pulsars can be described by five param- 
eters Wo, r], n, M, and N/,. Obviously, the three param- 
eters of the observed flux (the normalization, the index, 
and the cutoff) cannot constrain all eight parameters, 
but they can constrain some combinations of parame- 
ters. These constraints may be very useful if combined 
with results from other experiments, such as the obser- 
vations of protons, heavy nuclei, or diffuse gamma rays. 

Another concern is the reliability of constraints com- 
ing from the local e+e^ flux. Ideally, we would like to 
constrain the parameters in the models, but since the 
local distribution of pulsars is fundamentally random in 
nature, there is a possibility that we can only constrain 
the properties of particular pulsars without getting any 
information about the general population. The reason 
why we think our approach is sensible is the following. 
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At high energies, the flux from pulsars will depend sig- 
nificantly on the properties of individual pulsars (and we 
can use this region to prove that the observed flux is 
due to pulsars), but at low energies the flux is well ap- 
proximated by the continuous distribution flux and the 
properties of individual pulsars are relatively unimpor- 
tant. Thus, we propose using the observed spectrum at 
intermediate energies 100 GeV < E < 500 GeV as a test- 
ing ground to study the general (or averaged) properties 
of e~^e^ injection from pulsars. 

In the following we will fit the continuous distribu- 
tion flux derived in Eq. ([ST]) to the Fermi and PAMELA 
data simultaneously. In these fits we substitute the 
usual pulsar injection cutoff M ^ 10 TcV by the prop- 
agated (or statistical) cutoff Mgtat and treat it as a fit 
parameter. The other fit parameters are bo = (0.5 — 
3) X 10^^^ GcV^^ s^^ (these values are based on the 
energy densities of radiation and magnetic field within 
few kpc from the Earth), S ^ (0.3 - 0.6) [sB]. The in- 
jection index n w (1 — 2) and the conversion efficiency 
rjWo ^ — 10^^ erg are discussed in Appendix!^ In 
the fits we use Dq = 100 pc^kyr"^, iVb — 1.8 kyr^^ and 
Rgai — 20 kpc. The best-fit parameters are: 6 — 0.48, bo 
= 1.05 X 10-16 GeV-i s-\ r]Wo=A.9 x 10^* erg, n = 1.7, 
and Mstat = 1.0 TeV. The corresponding fluxes are shown 
in Fig. El 

In Fig. 10 we plot the 68%, 95% and 99% CL of these 
parameters. Every contour plot is obtained by varying 
two parameters while keeping the rest fixed at their best- 
fit value. We can see that for any value of the energy loss 
coefficient bo within the region chosen, there is a value of 
diffusion index S that can provide a fit within 95% CL. 
Additionally, a higher value of n suggests a lower value 
of 6, in agreement with the calculation of the e'^e~ flux 
from a continuous distribution of pulsars where we expect 
that F - £:-»-('5+i)/2 from Eq. Assuming a high 

value for the injection index (n > 1.6), the Fermi and 
PAMELA data could be fltted by a relatively large region 
of values of the propagation parameters (0.3 < S < 0.6 
and6o < 3 x lO^^^GeV^^s^i). On the other hand, values 
of n < 1.3 do not seem to give a very good fit to the data, 
with any combination of propagation parameters. Also if 
the total energy converted to e+e^ is smaller than 2x 10*^ 
erg, then for a pulsar birth rate of 1.8 per kyr regardless 
of the value of n the Fermi and PAMELA data cannot 
be explained by the continuous distribution of pulsars. 

To show the robustness of our procedure we applied the 
same analysis for two different backgrounds, a power law 



Pulsar 
Distribution 


bo 


5 


n 


V 


Oont. dist. Bi 


0.5-3.0 


U.30-0.D0 


1.35-1.95 


0.025-0.11 


Cont. dist. B2 


0.5-3.0 


0.30-0.60 


1.40-1.95 


0.020-0.085 


4kpc pulsars Bl 


0.5-2.4 


0.30-0.60 


1.25-1.80 


0.025-0.14 


4kpc pulsars B2 


0.5-2.0 


0.30-0.60 


1.15-1.95 


0.020-0.13 



TABLE I: Table of best-fit parameters within 95% CL. The 
values of bo are in units of lO-i*^ GeV-i s-\ In the first 
two lines, the ranges are given for a continuous distribution 
of pulsars. In the last two lines the equivalent ranges are 
given for the distribution of all pulsars within 4 kpc. Bl 
and B2 stand for the two backgrounds, simple power law and 
more realistic-conventional background [s^, respectively. We 
assume the mean initial rotational energy Wo = 10^'^ erg. 

background (Bl) and a more conventional background 
used in [5g (B2). In Table I, we present the 95% C.L. 
allowed region of values of the averaged ISM and averaged 
pulsar properties, using the two different backgrounds for 
the continuous pulsar distribution used. Alternatively, 
we used the properties of pulsars in the ATNF database 
with estimated distances d < A kpc. We present in Table 
I the derived constraints in the ISM averaged properties 
and universal pulsar properties n and rfWo- 

Better data on the flux of the high energy e+e^ and 
on the pulsar birth rate will make this analysis more suc- 
cessful in confining the parameter space that is relevant 
for the pulsar scenario. Tighter constraints of the back- 
grounds and the parameters of propagation through the 
ISM will be needed to confine the properties of pulsars 
themselves. 
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FIG. 10: The fits of the continuous distribution fiux to the ATIC and PAMELA data. The best-fit parameters are: S = 0.48, 

&o=1.05 X 10"^*^ GeV"^ tiWo=4.9 x 10*** erg, n = 1.7, and M^tat = 1.0 TcV. The contours show the confidence levels relative 
to the best-fit, = (x^ — X*)/<^o/, where x* is the chi squared for the best-fit paxameters. 
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